Coarse-graining a restricted solid-on-solid model 
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A procedure suggested by Vvedensky for obtaining continuum equations as the coarse-grained 
limit of discrete models is applied to the restricted solid-on-solid model with both adsorption and 
desorption. Using an expansion of the master equation, discrete Langevin equations are derived; 
these agree quantitatively with direct simulation of the model. From these, a continuum differential 
equation is derived, and the model is found to exhibit either Edwards- Wilkinson or Kardar-Parisi- 
Zhang exponents, as expected from symmetry arguments. The coefficients of the resulting continuum 
equation remain well-defined in the coarse-grained limit. 
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I. INTRODUCTION 

Driven, non-equilibrium interfaces have received much 
attention in recent years. Various models have been used 
to describe such systems. 0,0. These models may usu- 
ally be assigned to one of a small number of universal- 
ity classes, each characterised by a set of scaling expo- 
nents. To each universality class corresponds a contin- 
uum Langevin equation; such an equation may therefore 
be identified for each lattice model if the exponents are 
known, for example, from kinetic Monte Carlo (KMC) 
simulations. This procedure, however,faces difficulties 
when crossover effects are important 3]. To overcome 
this, Vvedensky Q has suggested usinga combination 
of an expansion of the master equation an d the 

dynamic renormalisation group (DRG) [H, Q to coarse- 
grain the resulting description, thereby directly obtaining 
a continuum Langevin equation in the large-scale, long- 
time limit. This program is a particular realization of a 
program suggested by Anderson [Io| . 

In this paper we implement this procedure for a re- 
stricted solid-on-solid (RSOS) model with both adsorp- 
tion and desorption. The master equation is expanded 
to obtain a set of discrete Langevin equations; these 
are then numerically integrated and compared to direct 
KMC simulations of the model, with which they are 
found to be in quantitative agreement. In the up/down 
symmetric case, an ad hoc procedure and the DRG both 
lead to the Edwards- Wilkinson (EW) equation 

—0(x,t) = is 2 — ^(x 1 t) + VDr ] (x,t) (1) 

as the macroscopic description of the model; here, rj is 
a zero-average Gaussian noise field with unit variance. 



For the asymmetric case, DRG arguments lead to the 
Kardar-Parisi-Zhang (KPZ) equation, 



—4>{x,t)=v 2 —4>{x,t) + \ 2 {—4>{ x ,t)\ +^D7 1 {x,t) 

(2) 

as the coarse-grained description. This is consistent with 
symmetry arguments and simulations pd|. The coeffi- 
cients of the coarse-grained continuum equations remain 
well-defined in the macroscopic limit. 



II. THE MODEL 

The model is a simple generalisation of the RSOS 
model introduced by Kim and Kosterlitz ^lj , and a spe- 
cial case of that introduced by Hinrichsen et al It 
is described (in one dimension for simplicity) by a vector 
h(t) of L time-dependent, integer- valued heights hi(t), 
i = 1,2,..., N; time is also a discrete variable. The 
dynamical evolution is given by the following rules: At 
each time step, an integer 1 < i < L is randomly cho- 
sen; with probability p the height hi is increased by 1, 
and with probability q attempt it is decreased by 1, pro- 
vided that the resulting configuration does not violate 
\h% — hi±i\ < 1. Finally, periodic boundary conditions are 
imposed. The allowed transitions, together with the as- 
sociated probabilities (equivalently, the rates) are shown 
in fig.n The transition rate from configuration h to h+r 
is given by 



{qS(n - l^A+^t-A-^) 
+ P 6(n + l)6{-A + hi)6(A-hi)\ 
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FIG. 1: Allowed transitions and associated transition rates. 



with Ti = U.jjti S ( r 3% 
± {hi±i - hi), and 

0(x) 



the discrete derivatives A ± h; 




(4) 



for integer (or zero) x. The temporal evolution of the 
probability density P(h, t) is given by the master equa- 
tion 



d_ 

dt 



P(h,t) 



d L r 



W(h-r;r)P(h-r,t) 
-W(h;r)P(h,t) 



(5) 



In Ref. |l3|, the special case q = 1, p = is stud- 
ied, and the KPZ equation derived as the coar se-g rained 
description. However, the treatment of Ref. [Tj| leads 
to ill-defined coefficients precisely in the coarse-grained 
limit. In particular, the authors find the coarse-grained 
description to be eqn. J2J) with V2 , A2 and D proportional 
to different powers of a parameter, a, which controls the 
degree of coarse-graining; indeed, their results are valid 
as a — > 0. Such a problem does not arise in our approach. 



III. A VARIANT OF THE VAN KAMPEN 
EXPANSION 

An often-used method for deriving Fokker-Planck 
equations for stochastic processes is the van Kampen ex- 
pansion [EE10- The method as described in ^,[3 is not 
directly applicable to systems such as the RSOS model 
because it requires a small parameter, l/f2, in which to 
expand. In effect, it assumes the existence of a macro- 
scopic law along with stochastic corrections to it, the rel- 
ative size of which is controlled by the expansion parame- 
ter 6]. In stochastic growth models, this is not the case; 
it is impossible to separate the time evolution into de- 
terministic and stochastic parts: the stochastic evolution 
is all there is. A more extended discussion of this point 
may be found in Refs. [Lll 115) . where it is shown that 
eqn. JSJ may nevertheless be approximated by a Fokker- 
Planck equation, which corresponds to the Ito Langevin 
equation 



dt 



K 



(i) 



(h) + v / 4 2) (h) % (t) (6) 



where summation over repeated indices is implied and 
the jump moments are defined by 



K 



{n) 



(h) 



dr r, 



,W(h;r). 



Eqn. © is essentially a set of simultaneous coupled 
Langevin equations. Obtaining a continuum version is 
conceptually similar to the (inverse of the) method of 
lines used to solve partial differential equations. How- 
ever, in the case of discrete interface models, the jump 
moments contain non-analytic step functions which 
necessitate a regularisation procedure. This is done in 
the next section. 



IV. DISCRETE LANGEVIN EQUATIONS 

To apply eqn. JSJ, the first two jump moments are 
needed. These can be calculated from W given by 
eqn. and are 

K^ih) = q8(A + hi)6(-A-hi) - pd(-A + hi)6(A-hi) 



and 



Kg\h)=[q6(A+h i )9(-A-h i 



where Sij is a Kroneckcr delta. 



(7a) 



■p0(-A + hi)0(A-hi)] 
(7b) 



Eqns. (|7ajl and (|7b|l are not completely specified by 
the lattice transition rules of sec. [H] but must be ex- 
tended to non- integer values of the hi. Furthermore, in 
deriving eqn. ©, the implicit assumption that H^(h;r) 
is analytic in h was made, which in turn implies that the 
extension of 6(x) to non-integers must be differentiable. 
This continuation of the to noninteger arguments 

is known as regularisation 0, 0, 0, . We choose the 
representation of 9 given by 



where, following 



0<x) = lim 6 A (x), 



exp {{x + a)/A) + 1 



In 



exp {x/A) + I 



(8) 



(9) 



The parameter a is the value of x below which 0(x) = 
and must therefore satisfy < a < 1 in order to agree 
with the lattice rules; it is otherwise free at this stage. 
A is effectively a "smoothing" parameter, and #a is an 
analytic function of x at the origin for all A > 0. More 
details may be found in 0. 

In some cases, such as the Wolf- Villain model analysed 
in [l^, the lattice rules can be used to infer the value of 
a. Here, simple arguments 01 show that a = leads to 
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L = 100 with q = 1 and p = 0; results from a KMC sim- 
ulation and numerical integration of the Langevin equa- 
tions for each of a = 0.6, a = 0.7 and a — 0.8 are pre- 
sented. It is evident that a w 0.7 gives a result that 
agrees closely with the KMC simulation. To test this 
agreement in the case q — p = 0.5, we simulate a system 
of size L = 1000; fig. shows the early-time behaviour. 
It again shows that a w 0.7 gives the best agreement be- 
tween KMC simulations and the stochastic formulation. 

In conclusion, the system of discrete Langevin equa- 
tions with the choice ~ 0.7 gives quantitatively the same 
results as the original model. 



FIG. 2: Variation of the evolution of the interfacial width with 
time for three values of a, as well as KMC simulation; q = 1 
and p = 0, L = 100. Both KMC and stochastic integration 
results averaged over 2000 runs. 
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V. THE CONTINUUM LIMIT 

The continuum (coarse-grained) limit of the symmetric 
case q = p may be taken using a simple ad hoc method, 
which is direct but physically opaque and cannot be ex- 
tended to the asymmetric case. The DRG is necessary 
to analyse the full model; it is considerably more gen- 
eral and conceptually clear, but correspondingly more 
complicated to carry out. The analysis presented in this 
section confirms the results of symmetry arguments and 
constitutes the final step in the identification of the SDE 
corresponding to the RSOS model. 



A. Ad hoc approach 

Using the regularisation of sec. II V I the step functions 
may be expanded about x = as 



FIG. 3: Early-time behaviour of the width for L = 1000 and 
q = p = 0.5 Both KMC and stochastic integration results 
averaged over 1000 runs. 



a discontinuous (non-analytic) dependence of (h) on 
h; this is against the spirit of the model, so a ^ 0. Un- 
fortunately, the argument cannot fix the exact value; this 
must be determined by comparing the results of numer- 
ical integration of eqn. JfjJ) to simulations of the lattice 
model. Accordingly, fig.[3shows a plot of w against t for 



Mx) = 5>„(A)x". (10) 

Replacing the set hi (t) by a function of a continuous argu- 
ment x, 4>{x, t) (which coincides with h%{b) for x = i), and 
expanding the discrete derivatives transforms the discrete 
set of SDEs of eqn. JSJ into a (partial) stochastic differ- 
ential equation for (f>(x, t). 

By using power counting arguments (see [LH for de- 
tails), we find that the relevant terms in this equation 
are, at most, 



dt 



v 2 



dx 2 



Vi 



dx i 



d 6 <j) <90<9 3 

OX° OX OX* \ OX 



x d , 

+ Ai3tt- -rf- 



dx V dx 



dx 2 



^2^4) (11) 



r 



where 77 is gaussian noise with unit variance and the co- q + p, while the rest are proportional to q — p; this is 
efficients may be found explicitly as functions of the b n . consistent with arguments based on up/down reflection 
The coefficients i>2, va, ^6, A13 and D are proportional to symmetry 
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In the symmetric case q = p, the only surviving 
terms are given by v 2 = \1v± = 360^6 = 6 6i, A13 = 
(36063 + 6i6 2 )/3 and D = 1. The coefficients b n may be 
found explicitly from eqn. @ . These satisfy the following 
inequality |l9j: 



MA) < 



A n (A) 



(12) 



where A n — > const as A — > + . 

A simple scaling approach would involve rescaling 
space, time and the field as x — > ex, t — > e z t and (j) —> e a <p, 
respectively, followed by taking the limit e — > + . In- 
stead, following 0, 0] , we will take the two limits e — > 
and A — > together. Putting A = e 5 with S > to be 
determined, the limit e — > + is well-defined only if 2 = 2, 
a = 1/2 and <5 = 1/2. The space-time scaling is diffusive 
(EW) , while the value of 5 has no direct physical signifi- 
cance. At the limit, the only surviving terms are the v-i 
term and the noise; all the others vanish. We therefore 
conclude that the coarse-grained limit of the symmetric 
RSOS model is the EW equation pol ]. 

Unfortunately, this direct coarse-graining procedure 
cannot be applied to the asymmetric case. In the next 
section, both the symmetric and asymmetric version of 
the model are studied using DRG arguments. 



B. Dynamic renormalisation group 

An equation corresponding to the limit of eqn. (|ll|l for 
q = p has previously been analysed by das Sarma and 
Kotlyar Jig; however, the DRG flow equation for A13 is 
not derived in Ref . [l]| ■ 

For this section, eqn. Ijlll) will be generalised to d di- 
mensions in the obvious way. Defining Kd = Sd/ {2tt) c1 , 



S d = 2n d ^/r(d/2) and g = \ 13 D/v%, we find, to first 
order in g (which plays the role of an effective coupling 
constant) 



d^ 2 
dAi 3 



= i/ 2 [z - 2 



1 d + 2 
2 SKd — 



6tf + 20 



d-2a) 



(13a) 
, (13b) 
(13c) 



(the non-renormalisation of D is a well-known conse- 
quence of the fact that the deterministic part of the equa- 
tion is conservative) . The flow equation for g is given by 



dg 
dl 



37 = 



K d 



2d(d + 2) 



(3d 2 + Ud + 28) 



(14) 



For any positive g (1 = 0) (which is the case here [T^). 
g = is an attractive fixed point |2lj |. This fixed point 
corresponds to z = 2 and a — (2 — d)/2, with A13 = 0, 
that is, it is the fixed point of the EW universality class, 
in agreement with the ad hoc approach of the previous 
section (as well as symmetry arguments and computer 
simulations). 

Application of the full machinery of the DRG to 
eqn. I|ll(l is unnecessary because definite conclusions may 
be reached by inspection of the equation. In terms of the 
Fourier transform of 4>(x, t), which we denote by the same 
symbol (j){k,uj), eqn. (|llf> (generalised to d dimensions) 
becomes 



G o (k,^)0(k) =»?(k,cj)+ f M 2 (q,k-q)0(q)0(k-q)+ f [M 3 (k, k - q - q', q, q'Mq)0(q')^(k - q - q')] ■ 

, 05) 



where the argument to has been suppressed for the 
bare response function is 

Go(k, uj) = f2k 2 — Vik* + vek e — iuj, 

and the two vertices are 

M 3 (ki, k 2 , k 3 , k 4 ) = A 13 [kj • k 4 ] [k 2 • k 3 ] (16) 

and 

M 2 (ki, k 2 ) = - A 2 (ki • k 2 ) + « 22 (k!) 2 (k 2 ) 2 
+ i Kl3 (k 1 .k 2 )[(k 1 ) 2 + (k 2 ) 2 ]. 



In the long-wavelength limit k — > the 2-vertex of 
eqn. (|17(l will be dominated by the A 2 (KPZ) term, with 
the Kij terms playing the role of higher-order corrections. 
Similarly, the 3-vertex is of higher order than the KPZ 
term; in addition, it has been shown to be irrelevant pre- 
viously. Therefore, the KPZ term A 2 determines the uni- 
versality class of the model. 

Since A 2 , «y cx (q — p), if q = p then M 2 = initially. 
If only vertices with an odd number of legs are present 
in the bare (unrenormalized) equation then vertices with 
an even number of vertices cannot be produced under 
renormalization. Therefore, if q = p the KPZ term is not 
present and the coarse-grained dynamics of the system is 
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described by the EW equation. If q ^ p, the dynamics is 
described by the KPZ equation. This result is consistent 
with symmetry arguments and simulation results [bH ]. 

VI. SUMMARY 

We have implemented a procedure suggested byVve- 
densky (and in more general terms by Anderson [l0|) to 
obtain macroscopic equations from microscopic models. 

Discrete Langevin equations are first derived and nu- 
merically integrated; they are found to be in quantita- 
tive agreement with KMC simulations of the underlying 
model. Next, these equations are expanded leading to 
a continuum Langevin equation, from which it is shown 
that the coarse-grained description of the model is the 
KPZ equation with the coefficient of the nonlinear term 
A2 vanishing in the symmetric case, so that the EW equa- 



tion is obtained. The coefficients appearing in the equa- 
tion are well-defined in the coarse-grained limit. 

The advantage of this procedure over the identifica- 
tion of the universality class by direct determination of 
the exponents from simulations is that slow convergence 
to the asymptotic regime is not a problem. In addi- 
tion, the DRG approach allows, in principle, investiga- 
tion of crossover effects (although this has not been pur- 
sued here) . Application of this procedure to other models 
would be a fruitful area for the future. 
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